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ABSTRACT 

Radiation  transport  modulates  the  spatial  frequencies  of  atmospheric  structures,  acting  as  a 
low  pass  filter,  which  causes  the  power  spectra  of  the  accumulated  radiance  to  have  different  power 
spectral  slopes  than  the  underlying  atmospheric  structure.  Additional  effects  arise  because  of  the 
non-stationarity  of  the  atmosphere.  The  SHARC  atmospheric  radiance  code  is  used  to  model  both 
equilibrium  and  non-equilibrium  radiance  and  radiance  fluctuation  statistics.  It  predicts  two 
dimensional  radiance  spatial  covariance  functions  and  power  spectral  densities,  PSDs.  Radiance 
power  spectral  slopes  for  paths  through  isotropic  Kolmogorov  turbulence  are  predicted  to  vary  from 
-5/3  to  -8/3  depending  on  the  length  of  the  path  through  the  turbulence.  The  input  gravity  wave  3-D 
covariances  and  PSDs  of  atmospheric  temperature  are  consistent  with  current  gravity  wave  theory, 
having  vertical  and  horizontal  power  spectral  indices  of  -3  and  -5/3,  respectively.  Altitude  profiles 
of  variances  and  correlation  lengths  account  for  the  non-stationary  of  the  gravity  wave  structure  in 
the  atmosphere.  The  radiance  covariance  and  PSD  power  spectral  slopes  differ  from  the  atmospheric 
gravity  wave  temperature  model  values  of  -3  and  -5/3.  These  modulations  depend  on  LOS 
orientations,  and  scale  lengths  of  the  sampled  altitudes  along  the  LOS. 


1.  INTRODUCTION 


This  paper  discusses  the  effects  of  radiation  transport  on  3-D  atmospheric  structures  induced  by 
gravity  waves  and  turbulence.  The  radiance  structure  of  atmospheric  backgrounds  is  induced  by 
local  temperature  and  density  fluctuations  in  the  atmosphere.  Radiation  transport  modulates  the 
spatial  frequencies  of  atmospheric  structures.  The  transport  acts  as  a  low  pass  filter,  and  causes  the 
power  spectra  of  the  accumulated  radiance  to  have  different  power  spectral  slopes  than  the 
underlying  atmospheric  structure.  Additional  effects  arise  because  of  the  non-stationarity  of  the 
atmosphere.  Horizontal  spatial  scales  of  the  atmosphere  structure  grow  with  altitude.  When  the 
structure  from  several  altitude  layers  contribute  to  the  radiance,  several  scale  lengths  are  present  in 
the  radiance  images.  Multiple  scales  can  introduce  breaks  in  the  radiance  power  spectral  slopes  and 
have  the  effect  of  lowering  the  spectral  slope  of  the  composite  power  law  PSD  at  spatial  frequencies 
between  the  high  and  low  frequency  structures.  The  Strategic  High  Altitude  Radiance  Code 
(SHARC)'  calculates  mean  line-of-sight  (LOS)  radiance  and  transmittance  values.  Its  structure 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

1998 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1998  to  00-00-1998 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Radiation  transport  effects  and  the  interpretation  of  infrared  images  of 

5b.  GRANT  NUMBER 

glilVllJ  miu  lUlUUlCllCC 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Spectral  Sciences  Inc, 4  Fourth  Avenue, Burlington,MA, 01803-3304 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distrihution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

14 

19a.  NAME  OF 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


2 

model  predicts  the  two-dimensional  spatial  covariance  function  of  the  radiance.  The  covariance 
function  and  its  Fourier  transform,  the  power  spectral  density  (PSD),  can  be  use  directly  in  sensor 
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models  or  in  image  synthesis  models  to  create  realizations  of  the  predicted  structure.  The  radiance 
statistics  and  images  are  non-stationary  and  are  explicitly  dependent  on  spectral  band  pass  and 
sensor  field-of-view  (FOV).  An  additional  quantity  derived  by  the  structure  model  is  a  LOS 
radiance  covariance  distribution  function.  The  radiance  covariance  distribution  function  acts  as  a 
weighting  function,  which  reveals  the  relative  importance  of  temperature  and  density  fluctuations  to 
radiance  variance  at  each  location  along  the  LOS.  The  radiance  covariance  distribution  function 
also  provides  the  distribution  of  scales  lengths  and  power  spectral  slopes.  The  SHARC  model 
predicts  that  radiance  power  spectral  slopes  depend  on  LOS  orientation  and  on  the  relative 
horizontal  to  vertical  scale  lengths  that  are  present  in  the  atmospheric  structures.  These  changes  in 
form  of  the  power  spectra  must  be  taken  into  account  when  making  inferences  about  the  atmospheric 
structure  from  the  analysis  of  infrared  imagery. 

There  is  considerable  evidence  for  isotropic  turbulence  in  the  troposphere"^’^,  and  for  gravity 
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wave  structures  in  the  mesosphere  ’  .  Recently,  evidence  for  thunderstorm  generated  gravity  wave 

o 

structures  in  infrared  images  in  the  upper  stratosphere  has  been  reported  .  The  SHARC  structure 
model  predicts  that  radiance  power  spectral  slopes  for  paths  through  isotropic  Kolmogorov 
turbulence^  are  modulated  to  vary  from  -5/3  to  -8/3  depending  on  the  length  of  the  path  through  the 
turbulence,  differing  from  the  isotropic  temperature  power  spectral  slope  of -5/3.  The  effect  is 
similar  to  the  wet  delay  interference  observed  in  the  radio  interferometric  measurements.  These 
experiments  measure  the  differential  phase  between  two  paths.  The  wet  delay  is  caused  by  water 
vapor  fluctuations  in  the  troposphere,  which  induce  refractive  index  fluctuations  along  the  LOS 
paths.  The  wet  delay  has  been  successfully  modeled"^’^  based  on  isotropic  Kolmogorov  turbulence. 
While  there  has  been  no  reported  modulating  effects  of  radiation  transport  through  gravity  wave 
structures,  the  emission  power  spectra  of  O2  airglow  and  of  the  OH  airglow  appear  to  be  steeper 
than  the  simultaneously  measured  rotational  temperature  power  spectra.  SHARC  contains  a  non¬ 
isotropic  gravity  wave  power  law  model.  The  predicted  radiance  PSD  power  spectral  slopes  differ 
from  the  power  law  temperature  and  density  fluctuation  models.  The  modulations  are  shown  to 
depend  on  LOS  orientation  and  on  the  relative  scale  lengths  along  the  LOS. 


2.  RADIANCE  STRUCTURE  CALCULATION  TECHNIQUE 

The  SHARC  model'  calculates  mean  LOS  radiance  and  transmittance  values.  Its  structure 
model  predicts  the  two-dimensional  spatial  covariance  function  of  the  radiance  and  a  radiance 
covariance  distribution  function.  The  covariance  function  and  its  Fourier  transform,  the  power 
spectral  density  (PSD),  can  be  use  directly  in  sensor  models  or  in  image  synthesis  models  to  create 
realizations  of  the  predicted  structure.  The  radiance  statistics  and  images  are  non-stationary  and  are 
explicitly  band  pass  and  sensor  FOV  dependent.  The  inputs  to  the  calculation  are  altitude  dependent 
temperature  and  density  profiles,  which  can  be  obtained  from  the  standard  atmosphere  generator 
(SAG)'*'  model  and  an  atlas  of  molecular  line  parameters  based  on  HITRAN-92."  The  non- 
stationary  altitude  dependence  of  the  temperature  fluctuations  are  modeled  by  profiles  of 
temperature  variance  and  horizontal  and  vertical  correlation  lengths.  Current  estimates  of  these 
profiles  are  taken  from  the  non-stationary  stochastic  (NSS)  model.  In  the  troposphere  and 
stratosphere  collisional  quenching  of  excited  vibrational  states  is  fast  compared  to  radiative  decay. 
Species  tend  to  be  in  local  thermodynamic  equilibrium  (LTE).  Above  the  stratosphere,  there  are  an 


insufficient  number  of  collisions  to  equilibrate  vibrational  energy  before  radiative  decay  occurs.  The 
vibrational  states  are  not  charaeterized  by  the  kinetie  temperature.  Under  these  eonditions,  the 
vibrational  states  are  in  non-local  thermodynamic  equilibrium  (NLTE).  SHARC  contains  NLTE 
models  for  both  mean  properties,  radianees  and  transmittances,  and  for  the  fluctuation  responses. 
Under  NLTE  eonditions,  the  response  of  population  ehanges  to  fluctuations  in  temperature  and 
density  are  determined  for  each  state  by  a  perturbation  model  which  uses  the  same  steady  state 
chemieal-kinetics  sehemes  that  are  used  for  the  quiescent  atmosphere.  Eluetuations  in  the  kinetie 
temperature  can  result  in  correlated  or  anti-correlated  responses  in  vibrational  state  population.  The 
magnitudes  of  the  responses  depend  on  the  loeal  ehemieal-kinetics  environment.  NLTE  emission 
depends  on  the  vibrational  temperatures  of  all  species  involved,  as  well  as,  the  kinetie  temperature. 

The  starting  point  for  determining  the  radianee  covariance  is  a  radianee  fluctuation  expressed  as  the 

2 

aceumulation  of  responses  to  local  atmospheric  fluctuations  in  temperature 


^AAp)  =  \drF^;,{r)M{r).  (1) 

includes  radiance  fluetuation  contributions  induced  by  the  rearrangement  of  population 
among  the  molecular  vibrational  states,  and  kinetic  temperature  fluctuations,  F°^{r) , 

F.,  (r)  =  ('’)«?■,.«  inienr)  +  (P) .  (2) 

j 

The  radiance  covariance  Cov^(j7,j7')  ,  where  p  and  j?' point  to  two  locations  in  the  sensor  image, 
can  be  expressed  as 


Cov^  {p,  p')  =  E[AL^^  ip'). 


(3) 


E[x]  is  the  expeetation  value  of  x.  Substituting  Equation  (1)  into  Equation  (3)  yields 


Cov^  {p,p')  =  j  dr^  drF^^  ir)F^;^  {r  ')Cov^  (r ,  r ') 


(4) 


It  is  assumed  that  the  radianee  eovarianee  and  the  3-D  atmospherie  temperature  eovarianee  are 
locally  stationary.  Under  this  assumption,  within  a  region  about  a  pixel,  p,  in  the  sensor  image,  the 

radiance  covariance  only  depends  on  the  lag,  Pp  =  p—p'  .  Likewise,  within  the  atmosphere  the 
temperature  covariance  is  only  a  function  of  the  lag,  Ar  =  f  -f' ,  about  the  point  f  . 

A  useful  intermediate  quantity  results  from  the  integration  over  the  LOS  r' .  The  result  is 
function  of  position  along  the  referenee  LOS,  and  deseribes  the  distribution  of  contributions  of  the 
atmospheric  structure  to  the  radiance  covariance.  The  resulting  radiance  covariance  distribution 
function  determines  the  LOS  spatial  distribution  of  the  radiance  covariance.  Eor  a  particular 


location  along  the  central  LOS,  the  integration  over  all  points  along  the  second  LOS  r '  is  equivalent 
to  integration  over  all  lags,  Ar  ,  equal  to  or  larger  than  the  minimum  lag  at  the  point,  Ar^ .  The 
minimum  lag  between  the  two  LOS  is  parallel  to  the  lag  veetor,  Aj? ,  in  the  image  plane,  as 
illustrated  in  Figure  1.  The  radiance  covariance  is  given  by 


Covl(p,Ap)  =  |drWL(r,Ary) 


(5) 


where  the  radianee  eovarianee  distribution  funetion,  WL(r,  Arjj)  ,  is  defined  as 

^ dr' .  (6) 

When  the  lag  in  the  image  is  zero,  the  minimum  lag  along  the  LOS,  Arjj ,  is  also  zero  .  Integration  of 
Equation  (5)  for  zero  lag  leads  to  a  LOS  variance  distribution  function,  or  weighting  function, 
IL^  (r  ,0)  .  The  radianee  variance  is  given  by 

<yl{P)  =  \drWi^{f,G).  (7) 


The  varianee  distribution  funetion  determines  how  eontributions  to  the  radianee  structure  are 
distributed  along  the  LOS  through  the  atmosphere.  The  integrations  over  the  two  LOS  through  the  3- 
D  atmosphere  cause  ehanges  in  the  form  of  the  radianee  covariance  funetion  relative  to  the 
temperature  covariance,  typically  yielding  a  different  spatial  frequency  dependence.  These  changes 
are  consequences  of  the  non-separability  and  non-stationarity  of  the  temperature  covariance  along 
the  LOS.  In  the  next  subsection  we  focus  on  the  first  integration  and  the  nature  of  the  resulting 
radiance  covariance  distribution  function.  In  the  subsequent  seetion  we  illustrate  the  effects  of  non- 
stationarity  in  the  proeess  of  evaluating  the  radianee  eovarianee  by  integrating  the  radianee 
eovarianee  distribution  funetion. 


3.  EVALUATION  OF  THE  RADIANCE  COVARIANCE  DISTRIBUTION 

FUNCTION 

In  this  seetion,  we  eonsider  several  limiting  eases  of  the  radianee  eovarianee  distribution 
function.  First,  a  special  case  of  a  single  narrow  layer  of  structure  is  considered.  Second,  an 
isotropic  temperature  covariance,  which  is  relevant  to  the  isotropic  turbulence  situation,  is 
considered.  Finally,  the  case  of  non-isotropic  structures  is  considered,  using  a  temperature 
covariance  model  that  describes  gravity  wave  structure  in  the  atmosphere. 


3.1  Very  narrow  altitude  layers  of  structure 

If  the  structure  is  dominated  by  a  single  narrow  altitude  layer,  the  integrand  in  Equation  5  may 
remain  approximately  constant.  If  the  radiators  in  the  band  pass  are  present  only  in  a  narrow 
altitude  layer,  then  (r ')  is  zero  except  for  a  short  segment  of  the  EOS,  A^ ,  located  at  the  center 
of  the  layer,  Vq.  The  integral  in  Equation  5  can  then  be  approximated  by 


(''o ,  )  ~  (^0  )Cov^  (ro ,  Arjj )  (8) 

In  this  approximation,  the  radiance  distribution  function  is  non-zero  only  at  the  same  altitude  layer 
as  the  temperature  fluctuation  and  the  radiance  covariance  is  given  by 

Cov^{p,Ap)^  As W, ip ,  Afjj )  (9) 

The  resulting  radiance  covariance  is  equal  to  the  atmospheric  temperature  covariance  in  the  direction 
parallel  to  the  image  plane.  The  conditions  for  this  to  be  valid  are  rather  stringent.  The  path  of  the 
LOS  through  the  layer,  A^ ,  must  be  short  compared  to  the  correlation  length  of  the  temperature 
structure  along  the  LOS,  in  order  for  the  temperature  covariance  to  be  constant  in  the  LOS  direction. 

This  will  become  more  apparent  when  paths  through  isotropic  and  general  non-separable 
structures  are  analyzed. 


3.2  Isotropic  temperature  structure 

The  second  special  case  is  that  of  an  isotropic  temperature  covariance.  This  is  important  for 
describing  isotropic  turbulence  structures  in  the  low  altitude  and  boundary  layer.  To  simplify  the 
analysis  we  use  an  approximation  suggested  by  Lindquist,  Kwon,  and  Nagy.  If  the  radiance 
fluctuation  amplitude  is  a  slowly  varying  function  of  location  along  the  LOS,  the  product 
(f ')  can  be  replaced  by  the  square  of  the  radiance  fluctuation  amplitude,  (r  ) ,  at  r  and 
taken  outside  of  the  integral  in  Equation  5,  resulting  in 

(f,  Af^ )  «  (f )J dr'Cov^ (r ,  Af )  .  (10) 

This  is  just  the  opposite  of  the  narrow  layer  approximation  made  above.  Since  the  temperature 
covariance  is  isotropic,  its  dependence  on  the  lag  is  independent  of  direction  and  is  only  a  function 

of  the  magnitude  of  the  lag,  |Af  |  =  -y/Ax^  +  Ay^  +  Az^  .  If  we  let  the  image  plane  be  the  XY  plane, 

the  lag  ,  Arjj ,  whose  magnitude  is  | Arjj  \  =  Ap  =  -^Ax^  +  Ay^  ,is  also  in  the  XY  direction,  and  the 

LOS  integration  is  along  z.  The  covariance  is  not  separable  along  the  LOS  and  the  resulting 
covariance  distribution  is  modified  in  the  XY  plane.  As  an  example,  a  simple  exponential 

covariance  function  with  unit  variance  and  scale  length,  CoVj  =  ,  can  be  integrated  to  give 

W,{f,Ap)  =  =  Fl,{f)pKpp)  (11) 

where  is  a  modified  Bessel  function  of  the  second  kind  of  order  1 .  Integration  through  the 
isotropic  3-D  covariance  yields  an  isotropic  covariance  distribution  function,  but  it  is  not 
exponential.  The  integral  can  be  found  in  compendiums  of  integrals  of  Bessel  functions^"^.  The 
exponential  can  also  be  expressed  as  e  ^  •  The  general  result  of  the  integration  in 

Equation  10  is  to  produce  a  new  Bessel  function  whose  order  is  increased  by  1/2.  These  functions 
are  the  covariance  functions  of  the  power  law  PSD’s.  If  rj  is  the  order  of  the  Bessel  function,  then 

the  corresponding  power  law  PSD  has  the  form  (1  +  k^  which  in  the  high  frequency  limit 

=>  .  In  the  exponential  example,  the  power  law  spectral  index  changed  on  integration  from  -2 

for  the  exponential  to  -3  for  the  radiance  covariance  distribution  function.  In  general,  a  change  in 


order  of  Bessel  function  by  one  half  corresponds  to  a  change  in  spectral  index  of  minus  one.  If  one 
uses  a  simple  Kolmogorov  turbulence  model^  of  spectral  index  -5/3,  the  integration  leads  to  a 
radiance  covariance  distribution  function  with  a  spectral  index  of  -8/3.  This  is  also  the  basis  for  the 
theoretical  findings  in  the  wet  delay  modeling  of  radar  interferometry."^’^  A  caveat  is  that  the 
analytical  integrations  are  over  all  space.  To  determine  the  actual  length  of  the  path  that  is  required 
for  modulation,  the  integration  was  performed  numerically  for  a  grid  of  interval  widths  from  zero  to 
several  scale  lengths.  The  results  of  these  integrations  are  shown  in  Figure  2.  The  results  for  the 
simple  exponential  case  are  instructive  for  other  power  law  covariances.  The  covariance  distribution 
is  modified  almost  completely  for  all  lags  that  are  on  the  order  of  the  length  of  the  path  through  the 
structure.  After  the  path  length  is  longer  than  a  few  scale  lengths,  the  covariance  distribution  is 
essentially  of  the  next  highest  half  order  over  its  measurable  range.  This  confirms  that  the  only  time 
that  the  approximation  in  Section  3.1  is  valid  is  for  lags  and  temperature  covariance  scale  lengths 
that  are  large  compared  to  the  path  length  through  the  structure.  In  addition  to  changes  in  spectral 
slope,  the  integration  leads  to  a  change  in  correlation  length.  For  the  simple  exponential,  the  scale 
length  and  the  correlation  length  are  the  same.  However,  for  the  modified  Bessel  function 

covariances  of  order,/?,  the  correlation  length  is  equal  to  C,  =  Vn  r(;7  + 1  /  2)lT{r])  where  3  is  the 
gamma  function.  In  Figure  2b  the  correlation  length  increases  from  1  for  the  exponential  covariance 
to  a  correlation  length  of  B/2  as  the  path  through  the  layer  increases.  The  integral  in  Equation  5 
can  only  be  replaced  by  the  integrand  times  the  path  length  of  the  LOS  through  the  narrow  region  for 
lags  larger  than  the  path  length  through  the  structure  and  for  temperature  correlation  lengths  in  the 
LOS  direction  that  are  longer  than  the  path  through  the  layer. 

It  is  instructive  to  evaluate  the  integral  in  Equation  1 0  using  Lourier  techniques.  Expressing  the 
3-D  temperature  covariance  in  terms  of  its  3-D  PSD  yields 


j  dzCovx(r,  Ar)  = 

Integrating  first  over  dz  yields  a  Dirac  delta  function^^  d{k^)m  spatial  frequency  along  the  LOS 
direction,  z.  The  result  is 

Wi^{r,ls.p)^^dzCoVj{f,ls.r)  =  ^dk^e  ‘''’'‘^dkye  (^k^  +k^)  ■  (13) 

The  integration  through  the  3-D  temperature  covariance  along  a  LOS  yields  a  2-D  radiance 
covariance.  The  2-D  covariance  is  the  Lourier  transform  pair  of  the  3-D  PSD  with  zero  spatial 
frequency  in  the  LOS  direction. 

In  the  special  case  of  the  Komogorov  isotropic  turbulence  model ,  for  example,  the  1-D,  2-D 
and  3-D  power  law  spectral  indices  are  -5/3,  -8/3,  and  -1 1/3,  respectively.  The  2-D  temperature 
PSD  is  obtained  from  the  3-D  temperature  PSDby  integrating  over  all  of  the  spatial  frequencies  in 
the  third  dimension.  The  summing  over  all  spatial  frequencies  in  the  third  dimension  leads  to  the  2-D 
temperature  PSD  with  its  larger  spectral  index.  The  effects  of  the  LOS  integration  is  to  remove  all 
high  frequencies  in  the  third  dimension.  LOS  integration  generates  a  2-D  radiance  PSD  with  the  3-D 
temperature  spectral  indices.  The  radiance  covariance  distribution  function  for  the  Komogorov 
turbulence  has  a  2-D  spectral  index  of-1 1/3  and  a  1-D  index  of -8/3.  The  2-D  spatial  frequency 
distribution  observed  in  the  image  plane  is  the  3-D  spatial  frequency  distribution  of  the  atmosphere 
in  that  plane  for  zero  spatial  frequencies  along  the  LOS  The  LOS  integration  eliminates  the  3-D 


j dzj dk^e”'*"^^  jdkyC  ''"^^PSDj.dC 


kz  +kx  +ky  ) 


(12) 


spatial  structures  that  are  correlated  with  high  spatial  frequeneies  along  the  LOS.  What  is  retained  in 
the  image  plane  is  the  accumulation  of  the  high  and  low  spatial  frequency  structures,  which  are 
correlated  to  the  low  spatial  frequencies  in  the  LOS  direction  .  It  is  the  selective  retention  of  the 
low  spatial  frequency  structure  in  the  LOS  direction  that  causes  selectiveness  in  the  image  plane 
when  the  LOS  direetion  in  the  3-D  eovarianee  is  non-seperable  from  image  plane.  Equations  12 
and  13  above  show  that  the  low  pass  filtering  properties  of  radiation  transport  is  general  and  not 
restricted  to  power  law  models. 


3.3  Non-Isotropic  Gravity  wave  Model 

The  power  law  temperature  eovarianee  model  eontained  in  SHARC  is  eonsistent  with  eurrent 
gravity  wave  models.  At  any  point,  f  ,  along  the  LOS  the  temperature  eovarianee  is  assumed  to  be 
separable  into  horizontal,  p  ,and  vertieal  ,z,  eomponents 


CoVj  (r ,  Ar  )  =  Cov^  (r ,  Ap)CoVy  (r ,  Az)  . 


(14) 


The  vertical  component  has  a  spectral  index  of -3.  The2-D  horizontal  eovarianee  is  taken  to  be 
isotropic  with  a  one  dimensional  speetral  index  of -5/3.  There  are  two  speeial  situations  whieh 
arise  in  the  evaluation  of  the  radianee  eovarianee  distribution  which  can  be  analytically  determined. 
We  utilize  approximations  whieh  are  not  made  in  the  SHARC  model,  to  simplify  the  analysis  and 
provide  insight  into  the  results  generated  by  the  more  aecurate  model.  For  long  range  nadir  or  zenith 
viewing,  where  the  LOSs  are  approximately  parallel  and  vertical,  the  LOS  integrations  are  in  the 
vertical  direction  .  In  this  situation,  the  horizontal  covariance  does  not  depend  on  the  integration 
variable  ,  the  vertical  altitude  lag,  Az ,  and  Equation  5  can  be  approximated  for  nadir/zenith  viewing 


Ify  (r ,  Ar^  )  «  Ify  (z,  Ap)  =  Cov^  (z,  A/?)J  (i(Az)Ffy  (z)Cov^  (z,  Az)  =  Cov^  (z,  AyC»)(y(z)  .  (15) 


The  integration  over  the  altitude  lag  leads  to  an  altitude  weighting,  (y(z),  of  the  horizontal 
covariance.  This  is  a  different  result  from  the  isotropic  case,  since  the  radiance  variance  distribution 
has  the  same  power  speetral  slope  and  eorrelation  length  as  the  horizontal  temperature  eovarianee  at 
the  altitude  z.  Unlike  the  situation  for  isotropie  turbulenee,  nadir  views  through  the  separable  gravity 
wave  model  would  have  a  PSD  with  spectral  slope  close  to  -5/3.  While  the  same  bias  toward  low 
frequeney  strueture  in  the  EOS  direetion  exists  in  both  models,  by  being  separable  the  horizontal 
spatial  frequencies  are  unaltered  by  the  bias. 

Another  approximate  situation  oecurs  if  the  EOS  is  in  the  horizontal.  In  the  ease  of  horizontal 
viewing,  ignoring  earth  curvature,  the  EOS  path  is  through  the  horizontal  isotropie  strueture  and  the 
vertieal  temperature  eovarianee  ean  be  faetored  out  of  the  integral  resulting  in  the  radianee 
eovarianee  for  horizontal  viewing 


JV^(x,Ay,Az)  =  F^^(x)CoVy(x,Az)j ci(Ax)Cov^(x,Ap)  =  a)(x)CoVy(x,Az)CoVjj(x,Ay)  .  (16) 


Here  the  x  axis  is  taken  to  be  the  horizontal  direetion  for  the  LOS.  Sinee  the  temperature  eovarianee 
is  separable  ,  the  radianee  eovarianee  distribution  funetion  is  also  separable  with  a  vertieal 
eomponent  which  is  not  modulated  with  a  power  spectral  index  of  -3  and  an  unchanged  correlation 
length.  The  horizontal  component  is  modulated  by  the  LOS  integration  through  the  2-D  isotropic 
structure.  The  radiance  covariance  distribution  horizontal  power  spectral  slope  is  modified  from  - 
5/3  to  -8/3,  and  its  correlation  length  is  increased  over  the  atmospheric  temperature  values.  Also 
within  the  flat  earth  approximation.  The  radiance  covariance  distribution  vertical  power  spectral 
slope  would  equal  -3.  The  flat  earth  approximation  is  severe  and  the  actual  treatment  of  limb  views 
is  reserved  for  the  general  case. 

The  path  integration  effects  for  general  viewing  geometries  have  to  be  evaluated  numerically. 
The  integrations  are  performed  using  the  SHARC  structure  model.  Within  this  model,  LOS’ s  are 
diverging,  based  on  the  observer  and  structure  geometry.  The  model  assumes  a  spherical  earth. 
Furthermore,  the  fluctuation  amplitudes  are  included  in  the  integrations  so  that  no  assumptions  are 
needed  about  their  constancy,  that  is  SHARC  numerically  integrates  Equation  5.  Based  on  the  results 
of  the  special  cases,  we  can  anticipate  that  the  radiance  covariance  distribution  function  will  have  a 
different  functional  form  from  the  3-D  temperature  covariance,  and  that  the  changes  will  be  a 
function  of  the  LOS  angle  through  the  structure.  In  addition,  there  are  strong  effects  based  on  the 
ratio  of  horizontal  to  vertical  correlation  lengths.  The  radiance  covariance  distribution  function  was 
evaluated  for  a  35  km  thick  layer  of  structure  between  30  and  65  km  in  altitude  for  a  grid  of  ratios 
of  temperature  correlation  lengths.  The  sensor  was  placed  above  the  structure  and  its  look  angles 
changed  so  that  several  local  off  nadir  angles  could  be  sampled.  The  covariance  distribution 
function  was  evaluated  at  40km.  The  local  off  nadir  angle  is  the  intersection  angle  of  the  LOS  with 
the  structure  at  40km.  The  radiance  power  spectral  slopes  changes  for  horizontal  to  vertical 
temperature  correlation  length  ratios  of  1,  5,  and  20  are  reported  in  Figure  3.  For  all  correlation 
length  ratios  the  radiance  covariance  distribution  power  spectral  slopes  are  within  -1.66  to  -1.8  at 
nadir  viewing.  Here,  the  sensor  vertical  and  horizontal  are  parallel  to  the  earth.  When  the  viewing  is 
nearly  horizontal  the  sensor  vertical  is  perpendicular  to  the  earth  and  the  power  spectral  slopes  of 
radiance  covariance  distribution  function  are  close  to  the  -8/3  and  -3  for  the  horizontal  and  vertical, 
respectively,  as  predicted  by  the  approximate  model.  For  intermediate  local  off-nadir  angles  the 
behavior  is  a  function  of  the  correlation  length  ratios.  For  large  horizontal  to  vertical  ratios,  the 
sensor  vertical  spectral  slope  is  close  to  -4  for  a  range  of  almost  60°  between  30°  and  90°  off-nadir 
angles,  the  vertical  slope  rapidly  falls  off  toward  -5/3  as  the  sensor  vertical  points  more  and  more 
into  the  atmospheric  horizontal  direction.  The  horizontal  power  spectral  slope  is  very  close  to  -2 
over  almost  the  entire  range  of  off  nadir  angles,  except  for  the  limb  and  nadir.  When  the  ratio  for 
temperature  correlation  lengths  is  1,  the  sensor  vertical  spectral  slope  has  a  narrow  minimum  of 
about  -3.5  at  60  to  70  degrees  off  nadir  while  the  horizontal  spectral  slope  is  close  to  -8/3  from  60  to 
90  and  steeper  than  -2  for  almost  all  off-nadir  viewing  angles.  These  grids  represent  the  anticipated 
range  of  values  of  power  spectral  slopes  that  may  be  measured  that  is  consistent  with  the  gravity 
wave  power  law  model. 


4.  RADIANCE  COVARIANCE 


The  final  step  in  determining  the  radianee  eovarianee  is  the  seeond  LOS  integration,  ie  that  of 
integrating  the  radianee  eovarianee  distribution  along  the  LOS.  If  the  LOS  is  for  limb  viewing,  and 
the  radianee  eovarianee  distribution  funetion  is  dominated  by  a  single  altitude  layer  with  a  single  set 
of  eorrelation  lengths,  the  radianee  eovarianee  will  be  proportional  to  the  radianee  eovarianee 
distribution  funetion  at  that  altitude.  The  integration  of  the  radianee  eovarianee  distribution  funetion 
ean  dramatieally  alter  the  shape  of  the  radianee  eovarianee  funetion  when  different  seales  are 
involved.  The  non-stationarity  of  the  atmospherie  strueture  ean  lead  to  a  radianee  eovarianee 
distribution  funetion  that  has  widely  different  eorrelation  lengths  at  different  ranges  or  altitudes 
along  the  LOS.  If  there  are  signifieant  eontributions  from  altitudes  with  different  eorrelation  lengths, 
the  integrations  ean  lead  to  breaks  in  the  speetral  slope  of  the  radianee  eovarianee  and  to  shallow 
speetral  slopes  at  spatial  frequeneies  between  the  largest  and  smallest  seales.  An  example  of  the 
effeet  is  illustrated  in  Figure  4,  in  whieh  two  PSD’s  with  the  same  speetral  index  of  -8/3,  but  with 
eorrelation  lengths  of  4  km  and  40  km  are  added  together.  The  blend  has  a  break  at  the  roll  over 
point  of  the  shorter  eorrelation  length  PSD.  The  blend  has  a  shallow  speetral  slope  until  the  40km 
PSD  dominates  at  lower  frequeneies.  The  produetion  of  frequeney  regions  of  shallow  speetral  slope 
has  been  used  by  the  image  syntheses  eommunity  to  generate  pseudo  fraetal  seenes  from  simple 
auto  regressive  proeesses.  Short  range  viewing  of  strueture  of  eonstant  seale  lengths  ean  exhibit 
the  same  effeets  beeause  strueture  in  the  near  field  will  appear  larger  than  strueture  at  a  distanee.  The 
effeets  of  non-parallel  LOS’s  in  the  FOV  generates  eovarianee  distribution  funetions  with  apparent 
eorrelation  lengths  that  are  large  in  the  near  field  and  smaller  at  long  range  even  for  aetual 
temperature  strueture  whieh  is  stationary. 


4.1  Example  BTH  off  Nadir  FOV  for  CO^IvA 

The  Mideourse  Spaee  Experiment  (MSX)  has  eolleeted  a  large  quantity  of  IR  imagery  in  narrow 
bands  around  the  heart  of  the  4.3  :m  CO2  band.  The  data  ineludes  low  altitude  limb  views,  below 
the  horizon  views,  and  nadir  seenes.  These  data  should  provide  a  wealth  of  information  on  the 
atmospherie  struetures  present  between  20  km  and  60  km  in  altitude.  To  illustrate  the  predietions  of 
the  forward  model,  ealeulations  were  performed  for  the  narrow  4.22-4.36  :m  band  for  a  LOS  elose 
to  the  horizon,  sensor  off  nadir  angle  of  59  degrees.  The  eovarianee  distribution  funetion  was 
numerieally  ealeulated  and  profiles  of  its  varianee,  loeal  eorrelation  lengths  and  power  speetral 
slopes  were  determined.  These  profdes  are  reported  in  Figure  5.  The  horizontal  power  speetral  slope 
varies  from  about  -7/3  at  the  higher  altitudes  to  -8/3  at  20  km,  while  the  sensor  vertieal  slope  varies 
from  -4  at  the  higher  altitudes  to  -3.3  at  20  km.  This  refieets  the  behavior  deseribed  above  of  large 
versus  small  temperature  horizontal  to  vertieal  eorrelation  length  ratios.  The  radianee  eovarianee 
distribution  funetion  horizontal  and  vertieal  eorrelation  lengths  refieet  the  eontinuous  inerease  in 
seale  lengths  with  altitude  deseribed  by  the  quieseent  input  temperature  model  profiles.  The  varianee 
distribution  shows  the  altitude  weighting  of  the  eovarianee  distribution  funetions  eontribution  to  the 
radianee  eovarianee.  The  radianee  eovarianee  funetion  was  eomputed  and  its  properties  determined. 
The  relative  standard  deviation  (sigma/mean)  of  5%,  sensor  horizontal  and  vertieal  eorrelation 
lengths  of  1 1 .7  km  and  2  km,  respeetively  and  horizontal  and  vertieal  speetral  slopes  of  -2.2  and  - 
3.6,  respeetively,  were  determined.  The  horizontal  power  speetral  slope  of  -2.2  refieets  the  effeets  of 


accumulation  of  contributions  over  several  correlation  lengths.  It  is  smaller  than  any  single  local 
contribution  of  the  radianee  eovariance  distribution  function.  A  comparison  of  the  numerical 
covariance  versus  the  analytical  power  law  function  computed  over  a  range  of  50  km  shows  no 
breaks  or  any  significant  deviations  over  that  range.  This  is  the  range  of  continuous  correlation 
length  changes  in  the  quiescent  NSS  temperature  statistics  profiles. 


5.  CONCLUSIONS 

The  radianee  strueture  of  atmospherie  baekgrounds  is  indueed  by  local  temperature  and  density 
fluctuations  in  the  atmosphere.  Radiation  transport  effects  modulate  the  spatial  frequencies  of 
atmospherie  structures,  acting  as  a  low  pass  filter,  whieh  modifies  power  law  speetral  slopes.  The 
non-stationarity  of  the  atmosphere  eauses  structure  of  several  scales  to  be  in  the  FOV.  This  is 
expeeted  to  be  particularly  prevalent  whenever  several  altitude  layers  eontribute  strongly  to  the 
covariance  distribution.  Multiple  scales  can  introduce  breaks  in  the  power  law  spectra  and  lead  to 
shallower  speetral  slopes  of  the  eomposite  power  law  PSD  at  spatial  frequencies  between  the  high 
frequeney  and  low  frequeney  struetures  in  the  FOV.  The  importance  of  various  altitude  eontributions 
along  the  LOS  ean  be  assessed  by  the  varianee  distribution  funetion.  The  radianee  power  speetral 
slopes  depend  on  LOS  zenith  angles  and  on  the  relative  horizontal  to  vertieal  seale  lengths  that  are 
present  in  the  atmospheric  structures.  These  effects  must  be  taken  into  account  when  making 
inferences  about  the  atmospherie  structure  from  analysis  of  infrared  imagery. 
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8.  FIGURES 


Figure  1 .  Two  lines-of-sight  illustrating  the  calculation  of  the  radiance  covariance  distrihution 
function  and  the  radiance  covariance. 
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Figure  2.  Covariance  distrihution  functions  (left  figure)  as  a  function  of  normalized  lag  for 
several  layers  widths  and  correlation  lengths  (right  figure)  as  a  function  of  LOS  layer  width 
for  isotropic  structure. 
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Local  Off-Nadir  Angle 


Figure  3.  Angular  dependence  of  the  spectral  index  for  several  ratios  of  the  horizontal  to 
vertical  correlation  lengths  for  sensor  horizontal  (solid  line)  and  vertical  (dashed  line). 
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Figure  4.  Generation  of  shallow  power  spectral  density  slopes  and  breaks  from  multiple  scale 
structures. 
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Figure  5.  Altitude  dependence  of  the  covariance  distribution  properties  in  the  MSX  band  B1 
for  59  °  off-nadir  viewing  geometry. 


